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Abstract 

In a recent paper 1] using a conjecture it is shown how one can calculate the effect of a weak 
stationary gravitational field on vacuum energy in the context of Casimir effect in an external 
gravitational field treated in 1 + 3 formulation of spacetime decomposition.. In this article, em- 
ploying quantum field theory in curved spacetime, we explicitly calculate the effect of a weak static 
gravitational field on virtual massless scalar particles in a Casimir apparatus. It is shown that, as 
expected from the proposed conjecture, both the frequency and renormalized energy of the virtual 
scalar field are affected by the gravitational field through its index of refraction. This could be 
taken as a strong evidence in favour of the proposed conjecture. Generalizations to weak stationary 
spacetimes and virtual photons are also discussed. 
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I. INTRODUCTION 



Zero point energy is a concept that its reality has always been the subject of somewhat 
intensive debate. Recently in the light of its important role in the formulation of the so 
called cosmological constant problem it has been discussed by people with opposite opinions 
2|, y). In this same respect, one obvious question need to be answered is its response to a 
gravitational field. In doing so in jl| it is shown how combining the concept of a spacetime 
index of refraction (in 1 + 3 formulation of spacetime decomposition) and virtual photon 
propagation in a Casimir apparatus could lead to a general solution in a weak stationary 
gravitational field in a coordinate-free manner. The conjecture employed there is a simple 
one which assumes virtual photons are affected by the spacetime index of refraction in the 
same way that real photons do when propagate in a weak gravitational field. It is shown 
that starting with this conjecture and assigning a frequency change uj — > — to virtual photon 
frequencies, one ends up with the same equation for the gravitational force on Casimir energy 
as for any other mass/energy. In other words it is shown that the (finite) Casimir energy 
couples to all inertia like any other mass/energy. In what follows, using QFT in curved 
spacetime, as a first step we study Casimir effect for a massless scalar field in a weak static 
gravitational field to explicitly calculate both frequency and finite energy changes to show 
that it gives the same results conjectured in l| - This, in our opinion, is a strong evidence 
that the proposed conjecture is correct. One should note that scalar field being massless, 
from geometical point of view, moves on a null geodesic so it is expected that frequencies 
of the corresponding virtual massless particles to be affected in the same way that virtual 
photons do. On the other hand for the energy density of the scalar field, from the studies of 
QFT in curved spacetime and specially Casimir effect, it is well known j|3] that it should be 
one-half the electromagnetic case on account of the fact that electromagnetic field has twice 
as many modes. We will use this result to obtain the electromagnetic Casimir energy in a 
weak gravitational field from that of a scalar field. In the last section we will present some 
general arguments as to why it is expected that the same method should in principle work 
for stationary spacetimes. 
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II. FERMI METRIC AND WEAK STATIC GRAVITATIONAL FIELDS 



To solve the field equations for a massless scalar field in a weak gravitational field the 
first thing one should specify is the coordinates in which this is going to be done. On the 
other hand since the above mentioned conjecture is stated for weak gravitational fields, in 
this section we discuss in some detail our choice of coordinates for the weak field of a static 
gravitational field which is Fermi coordinates. In |3J it is shown that coupling to gravity 
of Casimir finite energy for a weak gravitational field depends on the coordinate system 
employed. Analyzing this difficulty it is shown that it basiclly arises from the simple fact that 
relations between coordinate increments and physical distances depend upon the distance 
from the gravitating center in most common coordinate system. Then it is discussed why 
the calculations based on Fermi coordinates are the ones that should be taken as physical. 
The main reason being the fact that the Fermi coordinate is the closest generalization to a 
weak gravitational field of an inertial coordinate system in flat space . So the metric of a 
weak static gravitational field, representing that of Earth or any other spherically symmetric 
mass, a distance z above its surface of radius R, is chosen to be the Fermi metric 

g 00 = -(l + 2gz) , fiv = <V > ^v=l,2,Z (1) 
with g = the gravitational acceleration. In brief this metric represents the proper 



2]. In what follows we 



reference frame of an accelerated observer in a curved spacetime 
will employ the basic quantum field theoretic approach also used in |9[ for the same problem 
but with the main difference that, in effect, the calculations are done in the weak field limit 
of Schwarzschild spacetime in isotropic coordinates. Of course by the above argument on 
the importance of the coordinates employed in such a calculation it is expected that our 
results to be at variance with those given in isotropic coordinates. We will compare the 
rersults and elaborate on this important fact in the last section. 

III. SOLUTION TO THE SCALAR FIELD EQUATION AND MODE FREQUEN- 
CIES 

Field equation for a massless scalar field $(2^) in a curved background outside the matter 
distribution is given by 

da[V = Bg ab d b $(a*)] = 0- t Q = detg ab ; a, b = 0,1, 2, 3. (2) 



To quantize a scalar field in a curved spacetime one defines the following generalization of a 
flat space scalar product |5| 

($i,$ 2 ) = -i I $i(x)0a$;(z)[-0E(x)]*ra°dE (3) 
7s 

in which n M is a unit timelike vector normal to the spacelike hypersurface £ (in our case z 
= constant hypersurfaces ) and rfS is the volume element in S. For a set of mode solutions 
of (EJ) orthonormal in the above product we have 

(® i (x),$ j (x)) = 6 ij 5(k i -k j ) (4) 

Now for a weak static gravitational field represented by the metric (CQ), the field equation 
(J2J) reduces to 

y°°^$ + V 2 $ + V(ZnV=fl).V$ = (5) 

Our main goal here is to solve the above equation for modes inside the casimir apparatus. 
For simplicity we consider the case in which the plates are orthogonal to the radius along 
the z-axis in a local flat spatial metric of ([I]) with its origin at the lower Casimir plate. Due 
to the static nature of the spacetime and the x — y symmetry inside the cavity one could 
choose the scalar field to be of the following general (positive frequency) form 

$( x ) = Ce- iwt e lkxX e ikyV Z(z) (6) 

in which C is the normalization constant found through the orthonormality condition (j4j) to 
be 

Substituting (EJ) in equation (jSJ) leads to the following equation for the unknown function 

Z{z), 

Z"(z)+P(z)Z'(z) + Q(z)Z(z) = (8) 
in which ' = d z , k\ = k\ + k 2 y and 

P(*) = (0.W=fl) ; Q(z) = -(kl + oo 2 g 00 ) (9) 
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Now employing the common practice of writing Z(z) = U(z)W(z) for the above linear, 2-nd 
order differential equation we end up with 



W = exp ( * f P{z)dz) = (io) 
i.l (-g)4 



U" + S(z)U = ; S(z) = g - ipa _ V (f 1) 

Now it is well known that if S' <C 2\S\2 one could use the WKB approximation to solve 
the equation for U. It is an easy task to show that in our case (the weak static field) this 
condition is satisfied using the fact that g — — (1 + 2gz) where 

S' = % u 2 -?,^f (13) 

and the condition for the WKB approximation reduces to g <C u or equivalently ^0- <C 1 
which is obviously true. Therefore from WKB approximation we have the following general 
solution 

U = -p^ ^ciexp(i J a/ S(z)dz) + c 2 exp(—i J a/ S(z)dz)^j (14) 
or appropriate for a Dirichlet boundary condition 

U = %tn(A(z) + <M (15) 
A(z) = [ y/S{z)dz (16) 



Applying the Dirichlet boundary condition to the above wavefunction on the plates i.e 
Z(0) = Z(£) = or equivalently U(0) = U{£) = we have 

A(0) + (f) = ln (17) 
A(£) + O = mir (18) 

leading to the relation 

ft i 

A(£) - A(0) = / ^S{z')dz' = nn (19) 
Jo 

in which l,m,n are all integers. Now taking into account that in our case goo = 1 + e with 
e = 2gz <C 2k±, from (fT2l) we have 

S l l 2 ^{-uj 2 g™-kl) l l 2 (20) 



which substituted in (|T9|) gives 

f\-u 2 g m -kl) l l 2 dz' = mx (21) 
Jo 

Working this integral out to the first order in g£ we find 

u = (1 + |> (22) 



with wo = A;]_ + (Jj-) the field frequency in the a 
relation is nothing but the proposed relation in 
in the cavity 



Dsence of the gravitational field. The above 



1] for the frequency of the virtual photons 



w = wo(foo) 3 = {l + 2g <z >)?u « (1 + g < z >)u (23) 

to the first order in gz. It should be noted that to compare the above two equations one 
has to take < z >= f°* dz — I in other words the midway between and £ (average of z) 
on the Z-axis is taken to be the point where the index of refraction is calculated. 
It should be noted that the above calculations are up to the first order in g£ (or gz) and 
this is so because the Fermi metric is the first order correction to the line element near an 
observer's worldline unaffected by the spacetime curvature. Spacetime curvature comes into 
play only at the second order leading to the Fermi normal coordinates in the absence of 
acceleration and rotation . As pointed out in the introduction it should be noted that 
one should do the above calculation for a photon field to account for the electromagnetic 
Casimir effect. But since our scalar field is massless it also propagates on a null geodesic so 
it is expected that the same frequency relation holds for virtual scalar particles as it does 
for virtual photons. 



IV. CASIMIR ENERGY IN FERMI METRIC 

Employing the setting introduced in the begining of the last section the vacuum energy 
density for a scalar field is given by 

e vac = n a n b < 0\T ab \0 > (24) 

which in the present case reduces to 

e vac = — S„ / d 2 k ± T 00 [<S> n (k), $;(fc)] (25) 
-0 J 
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where T a b[(p n (k), = Q ya $ ,b — \g a bg cd Q 4 is the bilinear expression for T ab [5j reducing 



2 

to the following expression for the T 00 component 



T 00 = \c 2 Z\z) ju, 2 - gooik\ + (f) 2 } (26) 

The cavity being in a curved background its energy density is position-dependent (in our 
case z-dependent) and so we consider the average Casimir energy density given by 

(e vac ) = jZ n J d 2 k ± I J J-Too[$ n (fc), K(k)]dz} (27) 

which is a consequence of the fact that both the Fermi metric (CQ) and Too are only z- 
dependent. We first calculate the integral inside the braces 

J ±T m [Q n (k),** n (k)]dz = J {-goo)- l \c 2 Z\z)[uj 2 - gookl- goo{^) 2 ]dz (28) 

by approximating the function Z(z) 

z{z) 00 (- d s\z)y/± s[ni J y^ dz + fa) ( 29 ) 



in which 

OJ 2 



^S{z)dz= /(— k 2 ± )*dztt I (u 2 -k 2 ± -2gu 2 z)*dz 

-goo 



j \az + b)^dz « Vbz (30) 



where b = u 2 — k\ and a = —2gu 2 . In the same way 



- 5 S(z))^ = bHl + [ 9 - + ^}z) (31) 



So to the first order in gz 

& Z f 1 

Z{z) « Z (l - [| + ^]z)sin(v / 62) and — « Vbcotg(Vbz) - -{g + ^) 

Now substituting from (1321) into (12 8j) and after some manipulation we have 

f I^L dz ^ l C 2 Zl [(l-2gz)(l-(g+^)z)sm 2 (Vbz)[u J 2 -goo(k 2 ± + (^ ) 2 }dz 
J -goo 2 J 2b ' - Z 



(32) 



1 -C 2 Z 2 u 2 £[l-(Ag+^ £ -} (33) 



On the other hand substituting for Z(z) in (J7J) from ( 1321) we have 



^2 = 1 ? (OA) 

2(2^uZ 2 £[l-(2g+f b )l] 1 J 
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So that (133]) reduces to 



7, '° dz»-£-r= (35) 



~9m 2(2^)2 

Substituting this back into the equation for the average energy density we have 

(tvac) = j^n J d 2 k ± | J -^T m dz \ w 2 (27r) 2 £ Sn / ^ 
in which by substituting for o> from the previous section we end up with 

(e vac ) = (l + 9 4) (t^v^ I d 2 k ± uj ) (37) 



2 ' V. 2 ( 27r ) 2 ^ 

and finally employing the Schwinger proper-time method and any of the well known reg- 
ularization methods to calculate the integral over the transverse momenta and the infinite 
sum in the above expression respectively, we get 

(e vac ) ren . = (1 + y)eL = (1 + 9 < * >)&c (38) 

In which e° ac = — 14 4Q£ 4 is the renormalized vacuum energy density for a scalar field in a 
Casimir apparatus in flat spacetime [5]. Although the electromagnetic case could be treated 
in the same way but as in the flat case it is expected that the electromagnetic value for 
the energy to be twice the above value for a scalar field on account of the fact that the 
electromagnetic field has twice as many modes i.e, 

(e V ac) EM = (1 + y)(24 c ) = (1 + y)("^) = i 1 + 9 < * >)(~^i)- ( 39 ) 

n 

Again this is the same result found in |l( on the basis of the conjecture discussed in the 
introduction. It is also in agreement with the result found in js] using the gravitational 
energy formula for a static, weak field spacetime represented by the Fermi metric. 

V. DISCUSSION 

Using basic concepts of quantum field theory in a curved backround it is shown how the 
frequency and renormalized finite energy of virtual scalar particles in a Casimir apparatus 
in a weak static gravitational field are modified with respect to their flat space values. 
The calculations are done in a Fermi coordinate system which was shown to be the closest 

n 

generalization to a weak gravitational field of an inertial frame in flat space [8|. As pointed 
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out previously, the same calculations are also done in the weak field limit of Schwarzschild 

n n 

spacetime in isotropic coordinates |9J. On the other hand in |3[, without employing quantum 
field theory in a curved background, it is shown that using the gravitational definition of the 
energy-momentum tensor for the electromagnetic Casimir stress tensor results in different 
values for the gravitational energy and force in different coordinate systems. In particular 
the values obtained for the gravitational energy and force in isotropic coordinates are shown 
to be twice those obtained in Fermi coordinates. Obviously there are no claculations for the 
photon frequencies in because of the non quantum field theoretric method used there, but 



surprisingly enough the frequency and energy modifications found in [9] (formulas (6.10) and 
(6.23)) for scalar particles in isotropic coordinates are twice those obtained in the present 



article 



13(. This is expected from the arguments about the relevance of the coordinate 



system used to represent the weak field limit in [3|. It should be noted that the value 
obtained for the finite Casimir energy in the present article was a direct consequence of 
the value found for the particle frequencies. Generalising our field theorertic approach to a 
stationary spacetime one could go back to the 1 + 3 formulation of stationary spacetimes in 
which the spacetime line element is written in the following form 

ds 2 = g 00 (dx° - A a dx a f - dl 2 (40) 

where A a = — ^ is the so called gravitomagnetic potential and 

dl 2 = ^ a pdx a dx p = (-g a p + ^^-)dx a dx p a, (5 = 1, 2, 3 (41) 

is the spatial length element in terms of the three-dimensional spatial metric j a p. In the 
weak field and slow orbit (v <C c) limit the spacetime metric could be written in the following 
form M, 



2$ 2A.v N 9 , 9 „ 9 

ds 2 « 1 + — )c 2 dt 2 - dl 2 42 

c c 



2* , a 

In which = goo and v a = In some sense in the metric form (1421) the cross 

term has been absorbed into the time-time component. In other words, in these limits, 
it looks as if one could think of a stationary spacetime as a static one with a more com- 
plicated goo component leading to the correct form for the spacetime index of refraction 
(n = -7!= rj 1 - $ + A.v; c = 1). 

Generalizing the results of the previous sections to the case of a weak stationary spacetime 



(say of a rotating spherical mass) the natural approach would be the application of the quan- 
tum field theoretic methods to the case of a Casimir apparatus in a spacetime represented 
by a Fermi coordinate system which includes rotation [6] 

ds 2 = (1 + 2a^)dt 2 + 2{A^)dt - S^dx^dx" (43) 

in which = e flu i3X u u 13 and u 13 is the observer's angular velocity relative to inertial-guidance 
gyroscopes. But from the above argument one could transform ( 143]) to the following equiv- 
alent form 

ds 2 = (1 + 2a li x> Jt + 2A^)dt 2 - S^dx^dx" (44) 

Comparing this metric with equation ([T]) it suggests that the previous results for a static 
spacetime will be repeated but now with the above modified goo component. Of course this 
should be compared with the result obtained from the direct application of quantum field 
theory in a background spacetime represented by metric (1431) [4j. 
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